Charge sensing and controllable tunnel coupling in a Si/SiGe double quantum dot 
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We report integrated charge sensing measurements on a Si/SiGe double quantum dot. The quantum dot is 
shown to be tunable from a single, large dot to a well-isolated double dot. Charge sensing measurements enable 
the extraction of the tunnel coupling t between the quantum dots as a function of the voltage on the top gates 
defining the device. Control of the voltage on a single such gate tunes the barrier separating the two dots. The 
measured tunnel coupling is an exponential function of the gate voltage. The ability to control t is an important 
step towards controlling spin qubits in silicon quantum dots. 



Following a proposal to use gated quantum dots to host 
spin qubits in semiconductors,^ a series of recent experiments 
using GaAs/AlGaAs heterostructures have demonstrated the 
measurement and manipulation of quantum dot semiconduc- 
tor spin qubits. ^'^'"^'^'^'^'^ At low temperatures, the T2 de- 
phasing time for semiconductor spin qubits is usually lim- 
ited by hyperfine coupling between the electron and nuclear 
spins.^'^^'^^ Because Si and Ge have abundant isotopes with 
zero nuclear spin, a number of proposals for spin qubits 
in semiconductors have focused on these materials r^*i^*i^ 

and gated quantum dots have now been fabricated in 

Si/SiGe,i6a7a8j9,20 s^.^os SyStemS,^^.^^.^3,24,25,26,27,28,29 Si 

nanowires,^^'^^'^^'^^ and patterned donor layers.^"^ 

One of the important features of quantum dots is that they 
can be easily coupled to one another in the form of a dou- 
ble quantum dot. The tunnel coupling t between the two dots 
in a double dot is an important quantity, because several dif- 
ferent qubit operations can be performed by controlling this 
tunnel rate as a function of time. For approaches in which the 
qubit is a single electron spin, control of t enables the per- 
formance of SWAP operations between two qubits.- In a two- 
electron triplet- singlet qubit, controlling t enables single-qubit 
rotations. And in an architecture in which each logical qubit 
is composed of three spins in three different quantum dots, 
control of t alone enables universal quantum computation.^^ 

Here we report measurements of a Si/SiGe double quantum 
dot with an integrated quantum point contact (QPC) used for 
charge sensing. The lateral confinement of the quantum dot 
is tuned elecrostatically using metal top gates, and the device 
characteristics allow comprehensive control of the inter-dot 
tunnel coupling. We demonstrate a smooth transition from a 
well-defined double dot with weak tunnel coupling, through 
a strongly coupled quantum dot molecule, to a single, large 
quantum dot. Charge sensing with the QPC is used to probe 
the transfer of charge from one dot to the other. The sharp- 
ness of this charge transition depends critically on the inter- 
dot tunnel coupling, and fits to this transition allow the tun- 
nel coupling to be extracted over the range from a weakly to 
a strongly coupled double dot. We measure a change in the 
tunnel coupling of more than a factor of 10, corresponding 
to more than a 100-fold modulation in the exchange coupling 
between the quantum dots. 

A scanning electron micrograph of the device design is 
shown in Fig. la. The two-dimensional electron system 
(2DES) is located approximately 75 nm below the surface of 
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FIG. 1 : (a) Scanning electron micrograph showing the top gate ge- 
ometry for the measured device. The transport current through the 
dot is indicated by hot and the charge sensor current by /qpc- (b) 
Coulomb diamonds measured in transport through the quantum dot 
tuned to be a single dot. (c) Charge sensing signal from the cou- 
pled quantum point contact showing the periodic steps in /qpc result- 
ing from charge transitions on the neighboring single quantum dot, 
shown schematically by the dashed oval superimposed on the SEM 
image, (d) Charge sensing signal from a weakly coupled double dot, 
showing two different charge sensing amplitudes. 



the hetero structure. The carrier density and mobility of the 
2DES determined from Hall and Shubnikov-de Haas measure- 
ments are 5 x 10^^ cm~^ and 44,000 cm^/Vs respectively. 
The four gates labeled T, L, M, and R are used to create the 
quantum dot confinement, while gate Q in combination with 
gate L form the electrostatically coupled QPC. This device 
design incorporates several features that appear to be impor- 
tant for the results reported here. First, in order to minimize 
the dot size, and similar to Ref. 7, this device does not have 
the "plunger" gates that are often found in double quantum 
dot designs. The result is that each side of the quantum dot is 
smaller than previously published Si/SiGe quantum dots that 
have implemented a similar gate geometry. ^^'^^'^^ Second, the 
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side-gates L and R are quite narrow (~ 40% the width of 
the side-gate in Ref . 20) in order to reduce the screening be- 
tween the dot and the QPC to maximize the sensor response to 
charge transitions. As we show below, this quantum dot can 
be tuned to be a single or a double quantum dot by adjusting 
the voltage on the middle gate Vm • 

In the single dot regime, with small |Vm|, symmetric 
Coulomb diamonds are observed in the transport current 
through the dot /dot as a function of the source-drain volt- 
age VsD and Vr, as shown in Fig. lb. The last measurable 
diamond gives a single-dot charging energy Ec = e'^/C = 
880/ieV and a capacitance between gate R and the dot of 
17 aF. The charging energy corresponds to a total dot capaci- 
tance C = 182 aF, and this gives a lever arm Cr/C — 0.092. 
After this last diamond, the tunnel barriers between the dot 
and the leads are too opaque to identify further charge transi- 
tions using transport current through the dot, and charge sens- 
ing must be used. 

Figure Ic and d show the use of a nearby QPC for non- 
invasive measurement of the charge states of the quantum 
dot(s). As Vr is made more negative, electrons are removed 
from the dot, and the resulting decrease in the electrostatic po- 
tential of the dot causes a sudden increase in /qpc . In addition, 
there is an overall linear decrease in /qpc due to the direct cou- 
pling between gate R and the QPC, resulting in the saw tooth 
pattern of /qpc vs. Vr shown in Fig. Ic. The uniform height 
and periodic steps in /qpc indicate the presence of a single 
quantum dot. When the dot is tuned into a weakly coupled 
double dot, the charge sensor reports two distinct amplitudes 
for transitions on the left and right dots, as seen in Fig. Id. 
The left dot is physically closer to the QPC, and it produces 
a larger step ^ 50 pA, while a one electron change on the 
right dot yields (5/ ~ 15 p A. These charge signals on the QPC 
correspond to ~ 7 x 10~^ jh and 2 x 10~^ jh for the 
left and right dot respectively, in good agreement with obser- 
vations in similar sy stems. 

Figure 2a-c show the transition from a single dot to a well- 
separated double dot, as indicated by the charge sensing sig- 
nal dlqpc/dVR as a function of Vr and Vl, taken at three 
increasingly negative values of the voltage on the middle gate 
Vm' The numerical differentiation of /qpc as a function of 
Vr highlights the charge transitions on the dot which ap- 
pear as black or white lines in these plots. In Fig. 2a, with 
Vm = —0.5 V, the charge transitions form a single set of par- 
allel lines, indicating that the system is composed of a single, 
large quantum dot. The charge transitions are nearly sym- 
metric with respect to Vr and Vl, implying that the single 
quantum dot is coupled roughly equally to gate R and gate L, 
as expected based on the symmetric lithographic gate design 
(Fig. la). By making the voltage on gate M more negative 
(Vm = —0.7 V), the single dot is split into a strongly coupled 
double dot, where the left side is coupled dominantly to gate 
L and the right side to gate R, as shown by the bending of the 
transition lines in Fig. 2b. Further decreasing Vm to —1.0 V 
results in the weak-coupling regime and a well-isolated pair 
of quantum dots (Fig. 2c). 

The charge transition lines in Fig. 2c form the expected 
honeycomb pattern where each cell of the honeycomb cor- 
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FIG. 2: (a)-(c) dlqpc/dVR measured as a function of Vr and Vl 
for increasing \Vm\ (-0.5 V, -0.7 V and -1.0 V) (gray-scale shown 
in (c)). Three regimes of the quantum dot are evident from the dark 
transition lines: a single, large dot, a strongly coupled quantum dot 
molecule, and a weakly coupled double dot. (d) /qpc at a single 
vertex of the honeycomb diagram with a background plane removed, 
showing four charge states of the double dot. 



responds to well-defined electron occupations of the left and 
right dots.^^ Previous observation of double dot behavior in 
Si/SiGe was made in a structure that was designed to be a 
single quantum dot, and numerical modeling was required to 
understand how double dot behavior could arise in such a 
structure. In contrast, the structure used here and the data 
shown in Fig. 2 lead to a natural identif cation of the dot loca- 
tions: from the slopes and the spacings of the charge transi- 
tion lines in Fig. 2c, we find that the left dot is coupled ~ 2.2 
times more strongly to gate L than gate R (the more horizontal 
transitions), while the right dot is coupled ~ 2.0 times more 
strongly to gate R than gate L (the more vertical transitions). 

Figure 2d shows /qpc at one of the vertices of the hon- 
eycomb stability diagram with a background plane removed. 
Four distinct values of /qpc are visible, corresponding to four 
different charge states on the double quantum dot. Moving 
an electron from the left dot to the right dot results in an in- 
crease in /qpc- As shown in Fig. 3a and b, this increase ap- 
pears as a white line in a plot of dlqpc/dVR. In Fig. 3a, with 
Vm = — 1.0 V, the transition appears as a sharp white line, in- 
dicating a sharp interdot transition and weak tunnel coupling 
between the dots. In Fig. 3b, Vm = —0.825 V and the transi- 
tion appears as a broad and less intense white line, indicating 
strong tunnel coupling and increased charge delocalization be- 
tween the dots. 

Moving along the detuning diagonal in gate voltage space, 
indicated by the superimposed yellow line in Fig. 3a, causes 
a single electron to shift from one dot to the other. Measure- 
ments of this transition allow the extraction of a number of 
important quantities describing the charge sensor. The inset 
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FIG. 3: dlQPc/dVR measured as a function of Vr and Vl showing 
a single vertex of the charge stabihty diagram for the double quan- 
tum dot in the weak (a) and strong (b) coupling regimes (gray-scale 
shown in (a)). The superimposed yellow line in (a) indicates the de- 
tuning diagonal for the {N, M + 1) and (A^ + 1, M) charge states, 
(c) Inset: /qpc along the detuning diagonal of a weakly coupled dou- 
ble dot with the corresponding fit to the interdot transition given by 
Eq. IT]). Main plot: The interdot transition rescaled into units of ex- 
cess charge on the right dot at three different temperatures. The solid 
curves are fits using Eq. ^ with t = 0. Each of the data sets corre- 
sponds to an average of 15 sweeps. 



of Fig. 3c shows /qpc measured as Vl and Vr are simulta- 
neously swept to follow the detuning diagonal. On top of a 
linear background, due to the direct gate-QPC coupling, there 
is a step in /qpc when Vr ^ —0.2245 V due to the transition 
between the (TV, M + 1) and the (TV + 1, M) charge states. 

Following the method of DiCarlo et al.,^^ we treat the 
(A^, M + 1) and (A^ + 1, M) charge states as a two level sys- 
tem and fit the effect of the ground charge state transition on 
/qpc using: 

We)=/o + ^4tanh(^)+|e (1) 

Here, e is the detuning energy, Tg is the electron tempera- 
ture, a nd Q is th e ground and excited state energy splitting, 
Q = + 4t^, where t is the tunnel coupling, /q, SI, and 
dl/de are the background current, the sensor sensitivity, and 
the direct gate-QPC coupling, respectively. These three pa- 
rameters describe the system and are important for the extrac- 
tion of the tunnel coupling below. The detuning energy e is 
related to the gate voltage Vr by a lever arm a, e = aVR. As 
shown by the solid line in the inset of Fig. 3c, fits using Eq.[T] 
provides good agreement with the data. Using the results for 
the three fit parameters that characterize the QPC sensor, /qpc 
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FIG. 4: Excess charge on the right dot along the interdot transition 
taken with four different values of the double dot tunnel coupling t. 
Each data set is an average of 15 sweeps, and the solid lines are fits to 
the data using Eq. Inset: The extracted values of t as a function 
of Vm- Solid line is an exponential fit to the data. 



can be converted into units of excess charge on the right dot, 
(m) — M, and we report the data in the main panel of Fig. 3a 
and in Fig. 4 below in this way. 

When plotted versus gate voltage, the width of the interdot 
charge transition is dependent on three factors: Tg, a, and t. 
To determine Tg and a, we tune the double dot into the weak 
coupling regime, where t ^ 0, by making Vm increasingly 
negative until the transition is broadened only by temperature. 
In this limit, and with Vm fixed, raising Tg by raising the tem- 
perature of the dilution refrigerator broadens out the transi- 
tion, as shown in Fig. 3c. At low temperatures, the electron 
temperature is warmer than the refrigerator temperature, due 
to the balance between the cooling power supplied to the elec- 
tron system and heating from external sources. At higher tem- 
peratures, the electron temperature Tg approaches the refriger- 
ator temperature Tf. We fit the temperature dependence with 

a simple relationship between Tg and Tf,Te = ^^^T^^^T^, 

where Tq arises due to the effect of external heating. The pa- 
rameters To and a are determined by simultaneously fitting 
curves of the interdot transition in the weak coupling limit 
(t = and e = ft) at six different Tf using the formula, 

Im) -M =- I 1 + ^ tanh | ^ | | . (2) 

Three of the datasets and the corresponding fits are shown in 
Fig. 3c. We find a = 0.145 and To = 180 mK, correspond- 
ing to Tg = 182 mK at the base temperature of the dilution 
refrigerator. 

To determine t as a function of Vm, the extracted Tg and a 
are used to fit curves of the interdot charge transition. Fig. 4 
shows four such fits to the experimental data, with t as the sole 
free parameter. At small | Vm \ , the increased tunnel coupling 
results in a clear broadening of the interdot charge transition, 
corresponding to increasingly large t. 
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The largest uncertainty in the values for t extracted by this 
method arise from the uncertainty in the determination of a 
and Tg. As described above, a and Tg were determined from 
fits to the thermally broadened curves based on the assumption 
t = in the weak coupling limit. However, a good fit to those 
curves can be achieved with t as high as 6 /ieV, resulting in 
an uncertainty of 3% in a and 12% in Tg. The corresponding 
effect on the uncertainty in t is significant for the smallest t, 
while for large t the uncertainty is very small. 

We plot the resulting t as a function of Vm in the inset of 
Fig. 4. The data are well fit by an exponential curve. Changing 
Vm by ~ 200 mV changes t by more than a factor of 10, cor- 
responding to tuning the exchange interaction J by more than 
a factor of 100.^ In analogy with predictions of oscillations 
in the exchange coupling between donors, it is conceivable 
that t (and thus J) could oscillate as a function of Vm, due 
to the multiple valleys in Si. As shown in the inset to Fig. 4, 
within the experimental uncertainty we observe no indication 
of such an oscillation in the tunnel coupling t as a function 
of Vm . This smooth behavior could indicate either an effec- 
tively smooth quantum well interface in the active region of 
the device,"^ ^ as consistent with measurements of a large valley 
splitting in a QPC,^^ or it could be caused by valley relaxation 
in the presence of a rough interface."^^ 

Another important characteristic of quantum dots is their 
charge stability and noise behavior. In Fig. 2a-c, horizontal 
stripes of noise cut across the charging diagrams at several 
points. Based on their behavior as a function of the volt- 
ages on the gates defining the quantum dots and the QPC, 
we believe this noise is due to one or more charge traps be- 
neath gate Q. As is clear from Fig. 2, changing the voltage 
on gate L tunes the noise source in and out of resonance, and 
this behavior is consistent with previous reports that attributed 



some types of charge noise to the donor layer in semiconduc- 
tor heterostructures."^"^ Interestingly, this noise source seems 
to have little effect on the dots themselves. As can be seen in 
the figures, the lines corresponding to transitions in the quan- 
tum dot charge states pass through these noise stripes for the 
most part unchanged. The quantum dots themselves showed 
excellent long term stability, with charge drift less then 0.03e 
(which was the scan resolution) over a 7 day measurement pe- 
riod. This charge stability is consistent with other Si-based 
SETs which have shown drift of less than 0.0 le over a time 
period of weeks. 

In conclusion, we have demonstrated charge sensing mea- 
surements on a Si/SiGe quantum dot using an integrated quan- 
tum point contact. We have shown that the dot can be tuned 
into a single or double quantum dot. The device exhib- 
ited good charge stability, and the charge sensor measure- 
ments could resolve the interdot transition of a single elec- 
tron at fixed total charge. Broadening of the interdot transition 
caused by charge delocalization at reduced | Vm \ was used to 
extract t, and the tunnel coupling was found to have an ex- 
ponential dependence on this gate voltage. Control of t as a 
function of such a gate voltage is an important component in 
a number of approaches to quantum dot spin qubits. 

We acknowledge useful discussions about quantum dot gate 
geometries with L. M. K. Vandersypen and CM. Marcus. 
This work was supported in part by ARO and LPS (W91 INF- 
OS- 1-0482), by NSF (DMR-0805045), by DOD, and by DOE 
(DE-FG02-03ER46028). The views and conclusions con- 
tained in this document are those of the authors and should 
not be interpreted as representing the official policies, either 
expressly or implied, of the U.S. Government. This research 
utilized NSF- supported shared facilities at the University of 
Wisconsin-Madison. 



^ D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 (1998). 

^ J. Elzerman, R. Hanson, L. van Beveren, B. Witkamp, L. Vander- 
sypen, and L. Kouwenhoven, Nature 430, 431 (2004). 

^ J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, A. Yacoby, 
M. D. Lukin, C. M. Marcus, M. R Hanson, and A. C. Gossard, 
Science 309,2180 (2005). 

A. C. Johnson, J. R. Retta, J. M. Taylor, A. Yacoby, M. D. Lukin, 
C. M. Marcus, M. R Hanson, and A. Gossard, Nature 435, 925 
(2005). 

^ E H. L. Koppens, C. Buizert, K. J. Tielrooij, L T. Vink, K. C. 
Nowack, T. Meunier, L. R Kouwenhoven, and L. M. K. Vander- 
sypen, Nature 442, 766 (2006). 

^ E H. L. Koppens, K. C. Nowack, and L. M. K. Vandersypen, Rhys. 
Rev. Lett. 100, 236802 (2008). 

^ M. Rioro-Ladriere, T. Obata, Y Tokura, Y-S. Shin, T. Kubo, 
K. Yoshida, T. Taniyama, and S. Tarucha, Nat Rhys 4, 776 (2008). 

^ D. J. Reilly, J. M. Taylor, J. R. Retta, C. M. Marcus, M. R Hanson, 
and A. C. Gossard, Science 321, 817 (2008). 

^ A. Khaetskii, D. Loss, and L. Glazman, Rhys Rev B 67, 195329 
(2003). 

R. de Sousa and S. Das Sarma, Rhys Rev B 68, 115322 (2003). 
W. Yao, R.-B. Liu, and L. J. Sham, Rhys Rev B 74, 195301 (2006). 

B. Kane, Nature 393, 133 (1998). 



R. Vrijen, E. Yablonovitch, K. Wang, H. W. Jiang, A. Balandin, 
V Roychowdhury, T. Mor, and D. Divincenzo, Rhys. Rev. A 62, 
012306 (2000). 

M. Friesen, R Rugheimer, D. E. Savage, M. G. Lagally, D. W. 
van der Weide, R. Joynt, and M. A. Eriksson, Rhys Rev B 67, 
121301 (2003). 

A. Morello, C. C. Escott, H. Huebl, L. H. W. van Beveren, L. C. L. 
Hollenberg, D. N. Jamieson, A. S. Dzurak, and R. G. Clark, sub- 
mitted for publication, cond-mat/0904.1271v2. 
K. A. Slinker, K. L. M. Lewis, C. C. Haselby, S. Goswami, L. J. 
Klein, J. O. Chu, S. N. Coppersmith, R. Joynt, R. H. Blick, 
M. Friesen, et al., New J Rhys 7, 246 (2005). 
M. Sakr, H. Jiang, E. Yablonovitch, and E. Croke, Appl Rhys Lett 
87, 223104 (2005). 

T. Berer, D. Rachinger, G Rillwein, M. Muhlberger, H. Lichten- 
berger, G Brunthaler, and E Schaffler, Appl Rhys Lett 88, 1621 12 
(2006). 

L. J. Klein, D. E. Savage, and M. A. Eriksson, Appl Rhys Lett 90, 
033103 (2007). 

C. B. Simmons, M. Thalakulam, N. Shaji, L. J. Klein, H. Qin, 
R. H. Bhck, D. E. Savage, M. G Lagally, S. N. Coppersmith, and 
M. A. Eriksson, Appl Rhys Lett 91, 213103 (2007). 
E Simmel, D. A. Wharam, M. A. Kastner, and J. R Kotthaus, Rhys 



5 



Rev B 59, R10441 (1999). 

L. R Rokhinson, L. J. Guo, S. Y. Chou, and D. C. Tsui, Phys Rev 
B 63, 035321 (2001). 

A. Fujiwara, H. Inokawa, K. Yamazaki, H. Namatsu, Y Taka- 
hashi, N. M. Zimmerman, and S. B. Martin, Appl Phys Lett 88, 
053121 (2006). 

^4 G. M. Jones, B. H. Hu, C. H. Yang, M. J. Yang, R. Hajdaj, and 

G. Hehein, Appl Phys Lett 89, 073106 (2006). 

S.J. Angus, A. J. Ferguson, A. S. Dzurak, and R. G. Clark, Nano 
Lett?, 2051 (2007). 

N. M. Zimmerman, B. J. Simonds, A. Fujiwara, Y Ono, Y Taka- 
hashi, and H. Inokawa, Appl Phys Lett 90, 033507 (2007). 

H. W. Liu, T. Fujisawa, Y. Ono, H. Inokawa, A. Fujiwara, 
K. Takashina, and Y. Hirayama, Phys Rev B 77, 073310 (2008). 
G. R Lansbergen, R. Rahman, C. J. Wellard, I. Woo, J. Caro, 
N. Collaert, S. Biesemans, G. Klimeck, L. C. L. Hollenberg, and 
S. Rogge, Nat Phys 4, 656 (2008). 

E. R Nordberg, G. T. Eyck, H. Stalford, R. MuUer, R. Young, 
K. Eng, L. Tracy, K. Childs, J. Wendt, R. Grubbs, et al., in prepa- 
ration. 

Z. Zhong, Y. Fang, a Wei Lu, C. M. Lieber, and C. M. Lieber, 
Nano Letts, 1143 (2005). 

M. Hofheinz, X. Jehl, M. Sanquer, G. Molas, M. Vinet, and 
S. Deleonibus, Appl Phys Lett 89, 143504 (2006). 
Y Hu, H. O. H. Churchill, D. J. Reilly, J. Xiang, C. M. Lieber, 
and C. M. Marcus, Nat Nanotechnol 2, 622 (2007). 



F. A. Zwanenburg, C. E. W. M. V. Rijmenam, Y. Fang, C. M. 
Lieber, and L. R Kouwenhoven, Nano Lett 9, 1071 (2009). 

A. Fuhrer, M. Fuchsle, T. C. G. Reusch, B. Weber, and M. Y 
Simmons, Nano Lett 9, 707 (2009). 

J. Levy, Phys Rev Lett 89, 147902 (2002). 

D. P. Divincenzo, D. Bacon, J. Kempe, G. Burkard, and K. B. 

Whaley, Nature 408, 339 (2000). 

L. Dicarlo, H. J. Lynch, A. C. Johnson, L. I. Childress, K. Crock- 
ett, and C. M. Marcus, Phys. Rev. Lett. 92, 226801 (2004). 
W. G. Van Der Wiel, T. Fujisawa, S. Tarucha, and L. P. Kouwen- 
hoven, Rev Mod Phys 75, 1 (2002). 

N. Shaji, C. B. Simmons, M. Thalakulam, L. J. Klein, H. Qin, 

H. Luo, D. E. Savage, M. G. Lagally, A. J. Rimberg, R. Joynt, 
et al., Nat Phys 4, 540 (2008). 

B. Koiller, X. Hu, and S. Das Sarma, Phys. Rev. Lett. 88, 027903 
(2001). 

M. Friesen, M. A. Eriksson, and S. N. Coppersmith, Appl Phys 
Lett 89, 202106 (2006). 

S. Goswami, K. A. Slinker, M. Friesen, L. M. McGuire, J. L. Tru- 
itt, C. Tahan, L. J. Klein, J. O. Chu, P. M. Mooney, D. W van der 
Weide, et al., Nat Phys 3, 41 (2007). 
M. Friesen and S. N. Coppersmith, unpublished. 

C. Buizert, F. H. L. Koppens, M. Pioro-Ladriere, H.-P. Tranitz, 

I. T. Vink, S. Tarucha, W. Wegscheider, and L. M. K. Vander- 
sypen, Phys. Rev. Lett. 101, 226603 (2008). 



